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ABSTRACT 

Context. Transiting planetary systems allow us to extract geometrical information, e.g., the angle e between the orbital angular 
momentum and the stellar spin, that can be used to discriminate among different formation and evolutionary scenarios. This angle 
is constrained by means of the Rossiter-McLaughlin effect observed on radial velocity and can be subject to large uncertainties, 
especially for hot stars (T^ff > 6250 K). It is thus interesting to have an alternative method to constrain the value of the obliquity e and 
to detect companions that might have disturbed the orbit of the planet. 

Aims. We show how the long-term variations in the transit duration (TDV) can be used to constrain the obliquity of the stellar rotation 
axis. Our calculations may also be used to put an upper limit on the contribution of geometrical effects to the TDVs, thus allowing us 
to indirectly infer the presence of additional companions. 

Methods. We introduce a simple theory to describe the secular variations in the orbital elements and their effects on the TDVs with 
a general formulation valid for both oblique and eccentric systems. Parameters or orbital elements that cannot be directly measured, 
such as the longitude of the ascending node of the orbit, are avoided thus allowing us to perform a straightforward application. 
Results. We compute the expected TDVs for the presently known transiting systems, adopting their parameters found in the literature. 
Considering the capabilities of the present or next generation space-borne telescopes, we point out the systems that could be readily 
observed and discuss the constraints derivable on their fundamental parameters. 

Conclusions. Measured TDVs can be used to constrain the obliquity of the stars (and possibly of the planets in systems younger than 
10 - 100 Myr), giving information about the formation scenarios, the strength of the tidal coupling, and the internal structure of both 
the stars and the planets. Moreover, they can provide an indirect indication of other bodies, even with a mass comparable with that of 
the Earth, because they give rise to additional contributions to the nodal precession. 

Key words. Planetary systems - stars: rotation - planet-star interactions 



1. Introduction 

Planetary systems with transiting planets are a fundamental 
source of information on stellar and planetary properties, in par- 
ticular they allow us to measure the mass and the radius of the 
planets by removing the inclination degeneracy that affects sys- 
tems without transits. Moreover, the angle e between the orbital 
angular momentum and the stellar spin axis can be measured in 
such systems and enable us to discriminate among different for- 
mation and evolutionary scenarios. In the case of hot Jupiters, 
if the planet underwent inward migration by angular momentum 

exchange with a protoplanetary disc ( Lin et al 1996) . the stellar 

spin should be almost aligned with the orbital angular momen- 
tum, whereas if the small semi-major axis of its orbit were due 
to the Kozai mech anism (Kozai 1962) or planet-planet scatter- 
ing jRasio & Fordlll996l) . high obliquities would generally be 
expected. If the planet transits across the stellar disc, the radial 
velocity anomaly detected during th e transit, i.e., the so-called 
Rossiter-McLaughlin effect (see e.g. Triaud et al.ll20 1C/, and ref- 
erences therein) can be used to measure the projected obliquity 
A, i.e., the angle between the projections on the plane of the sky 
of the orbital angular momentum and the stellar spin. From the 
value of A, it is possible to derive the obliquity e if the inclina- 
tion is of the stellar spin to the line of sight is known (cf. Eq. |6] 
of Sect. Oi- 
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Since it has been empirically suggested that hot stars (7^,^ > 
6250 K) with hot Jupiters have high obliquities (IWinn et akl 
201Ci), it is very interesting to apply different methods both to 
constrain the obliquity of transiting systems in general and to 
look for the presence of additional companions that may have 
perturbed the orbit of the planet. 

In this work, we illustrate how the long-term variation in the 
transit duration (hereafter TDV) can be used to constrain the ob- 
liquity of the stellar rotation axis and to derive indirect evidence 
of other bodies in a planetary system. We study the general case 
of a system with both an oblique and eccentric orbit giving an ex- 
pression for the TDV that contains only measurable quantities. 
This is useful for a straightforward application of our results to 
the known transiting systems and to apply the method to ground- 
based and space-borne transit surveying and monitoring. 



2. Apsidal and nodal precessions 

Perturbations of the Keplerian orbit occur whenever the force 
field controlling the motion of the planet is not spherically sym- 
metric or doe s not follow the inverse -square dependence on the 
orbital radius. iMiralda-Escudel (l2002h considered two cases: the 
precession of the orbital plane of a circular orbit in a misaligned 
system and the precession of the periastron when the orbit is ec- 
centric, in the limit of a small eccentricity and assuming that the 
orbital plane coincides with the stellar equator In both cases, the 
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cause of the precession is the gravitational quadrupole moment 
of the star. The latter case was further extended to an arbitrar- 
ily large eccentricity by Pal & Kocsis (2008), independently of 
the physical mechanism causing the precession of the periastron. 
Iloriol (1201 ll) then proposed a set of uniform self-consistent ana- 
lytic expressions accounting simultaneously for the effects of all 
the perturbing accelerations on the duration and the mid-time of 
the transit. 

In the present work, we identify the known systems for 
which a TDV can be easily observed over a time interval of 
the order of a decade, so that they can be readily used to con- 
strain the spin-orbit alignment. Moreover, the apsidal precession 
in systems with an eccentric orbit may be used to estimate the 
Love number of the plane t that is related to i t s inte rnal density 
stratification, as shown by'Rago zzine & Wol3 ( l2009l) . 

We introduc e a sim ple approach, using the formulations of 
lMardling&Liiil(l200l for the secular variations in the orbital 
elements, merging the results of iMiralda-Escudel (|2002|) and 
iPal & KocsisI (l2008l) . to assess their effects on TDVs. In our for- 
mulation, parameters or orbital elements that cannot be directly 
measured, such as the longitude of the ascending node of the 
orbit, are avoided, allowing us a direct application of our theory. 

2.1. Rate of change in the orbital elements 

We assume a cartesian orthogonal reference frame :Ro with ba- 
sis vectors (i,j,k), the origin at the barycentre of the host star, 
and the k vector along the line of sight, so that the (i, j ) plane 
coincides with the plane of the sky. Following Mardling & LinI 
(1200 2). the orbital elements are expressed through the specific 
relative angular momentum vector h = r x f , where r is the po- 
sition vector of the barycentre of the planet, and the Runge-Lenz 
vector e, a vector in the direction of the periastron with magni- 
tude equal to the eccentricity of the orbit. By definition, the vec- 
tors h and e are always perpendicular, although their components 
in the adopted reference frame will change in time owing to the 
perturbing potential. The secular evolution in the elements is ob- 
tained by time-averaging over the orbit, using the true anomaly 
as the integration variable. The eccentricity e, the semimajor axis 
a, the inclination /, the argument of periastron co, and the longi- 
tude of the ascending node Q may be obtained from h and e 
through equations (21)-(27) of iMardUng & LinI (|2002|) . 

For a two-body system, the average rate of change in h and 
e can be written as the sum of the different perturbation contri- 
butions averaged over one orbit, i.e. 

(h) = <dQD,) + <dQD,> + (dxF,) + (dxFp), (1) 
<e) = <gQD,) + (gQDp) + <gTF,> + <gTFp> + <gGR>, (2) 

where the different subscripts denote the perturbations as fol- 
lows: 

- QDi is the effect of the quadrupolar distortion of the star, 
for / = s, and of the planet, for / - p, respectively; the 
quadrupole moment of each body arises from its spin cen- 
trifugal distortion as well as from the tidal distortion caused 
by the other body; 

- TFi is the effect of the tidal damping within the star, for / = s, 
and within the planet, for / = p, respectively; 

- GR is the effect of the (gravitoelectric) post-Newtonian po- 
tentiaQ 

3fl n e 

<gGR> = — 5—- rq, (3) 



' We neglect the usually much smaller grav itomagnetic effect, i.e., 
the Lense-Thirring precession (see Iloridl20l"ll) . 



where c is the speed of light, q = h x e is the unit vector 
orthogonal to the unit vectors in the directions of h and e, 
respectively, and n is the mean motion 



G(M, + Mp), 



(4) 



G being the gravitation constant, and and Mp the mass of 
the star and the planet, respectively. 

We focus on the evolution of the orbital elements that affect the 
observables of the transit on a relatively short timescale, i.e., of 
the order of a decade. Consequently, the terms dxFs, dxFp, gxF,, 
and gxFp related to the tidal effects are neglected, as they induce 
variations with typical timescales at least four orders of mag- 
nitude longefl Therefore, the angles of the stellar and planetary 
spins to the total angular momentum will be considered constant, 
as well as the values of the orbital eccentricity and semi-major 
axis. 

Considering that the planetary spin evolves to a coplanar 
state of pseudo-synchronization with the orbit within a typi- 
cal timesc ale of t ^ 10'^ yr fo r a planetary tidal quality factor 
Q'^- 10^ dLeconte et al. iBOTol) . we assume that the planet is al- 
ready in this state and thus the planet's angular velocity can be 
written as Wp = « h. The angular velocity of the star Ws can be 
described in by the angles is and (p^, which are respectively 
the inclination to the line of sight and the azimuthal angle, so 
that the corresponding unit vector can be written 



sin is cos ( 
sin /s sin ^ 
cos /c 



(5) 



The angular velocity of the star Ws is related to its equatorial ro- 
tational velocity v (assuming rigid rotation) by v - WsRs, where 
/?s is the radius of the star. The spin of the star Ls is defined as 
Ls = CsWs, where Cs is the tensor of inertia of the star Its trace 

Cs - aMfiR^ is the moment of inertia of the star, where a is the 
fractional gyration radius. 

We define e to be the angle formed by the orbital angular 
momentum vector and the spin of the star, viz. cos e = Ws ■ h, 
and let A be the angle formed by the projections of the orbital 
angular momentum and stellar spin on the plane of the sky. There 
is a simple relationship between 0s and A, i.e. 0s = + ^ + f ; 
moreover 



cos e = cos / cos + sin / sin cos A. 



(6) 



Applying thos e geometric relationships and inserting 
iMardling&Liiil 's Eqs. (48) and (49) into Eqs. (42) and 
(43), and then using their Eqs. (29), (30), and (3 IX the 
quadrupolar distortions of the star and the planet and the general 
relativity effect cause an orbit-averaged variation in the orbital 
elements given by 



<5> 



-As V^s cos e sin is sin A, 

2 COS e(cos is - cos / cos e) 



= -AsWt 



sin ; 



= As 



, GMp 
-^(3cos2e-l) + ^/o(e) 
2 



(7) 
(8) 
(9) 



n + — ^/o(e) 



3fl « 



2„3 



1 



1 



- <— )cosi, 
at 



- Only in a few systems, e.g., CoRoT-11, might the tidal interaction 
induce variations in the orbital elements with observable consequences 
over a time interval of ~ 30 y r if the stellar tida l quality factor Q[ is 
smaller than (3 - 5) x 10* (see iLanza et alj201ll) . 
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where 



Me) = 15 



1 + + 



(1 



(10) 



and As and Ap are factors that depend on the perturbing acceler- 
ations caused by the stellar and planet bulges, respectively, i.e. 



n(l -e2)2 + M, 



MAIR, 



P n(l -e2)2 \^ ^ mJ[ a ' ' 



(11) 
(12) 



where and 7?p are the radius of the star and of the planet, 
respectively, is the apsidal motion constant of the star, and 

- 5^L, where is the Love number of the planet. For the hot 
Jupiters, t he Love number is expect ed to be in the range ~ 
0.1 - 0.6 jRagozzine & WoHll2009l) . so we assume that ki^ - 
0.3. The values of k^^ are computed as a function of the mass 
and the effe ctive temperature of the star using the stellar models 
bv lClaret1(li995b assuming a solar composition. For transiting 
hot- Jupiters, we usually find that Ap/As ^ 100, which implies 
that the planetary distortion is predominant for orbital periods 
shorter than ~ 3 days, allowing a measurement of the planet 
Love number (see Ragozzine & Wolf 2009). 

For polar and equatorial orbits (e - 90° and e - 0°, i.e. A - 
0°), there is no nodal precession because the torque applied to the 
planetary orbit by the stellar quadrupole vanishes. In contrast, 
the perturbation to the inverse-square law of the force introduced 
by the quadrupole term of both components and the relativistic 
effect cause an apsidal precession for all values of the obliquity 
e. 



2.2. Variation in the transit duration 

We denote the duration of the transit by H, which is defined as 
the interval between the times when the centre of the planetary 
disc intersects the limb of the star during the ingress and the 
egress, respectively. Neglecting the curvature of the projection 

of the orbit because of its inclination, H ^ where AB is the 
chord travelled by the centre of the planet during its transit on the 
stellar disc with a tangential velocity Vtan. We define the angle y 
so that AB = 2Rs cos y, thus 



2R, 

H cosy. 

^tan 



(13) 



The tangential velocity can be expressed as (iPal & Kocsi J2008h 

1-l-ecosv 

Vtan = an — , (14) 
Vl -e2 

where v is the true anomaly of the planet. We assume that the 
planet has a constant velocity during the transit equal to the ve- 
locity at mid-transit. The true anomaly v and the phase angle 6 
are related by - v + a> - n/l. When / = 90°, the primary 
minimum occurs when 6-0° and v = n/l - oj. For transiting 
systems, / - 90°, thus we can apply these relationships without 
introducing an appreciable error and obtain 



1 + e sin <jj 
Vtan = an — , and Vta 

Vl -e2 



e cos CO 

an to. 

Vl -e2 



(15) 



We note that for nearly circular orbits (i.e. when the eccentric- 
ity and/or the Lagrangian orbital elements ecosw and esinw 



are zero within a few cr), the true anomaly and a> are meaning- 
less, but the phase angle 9 has a well-defined value. In this case, 
the tangential velocity at the time of each transit do not change 
within the uncertainties. We define b' to be the projected dis- 
tance between the centre of the planet and the centre of the star 
at mid-transit, i.e, the impact par ametef, b' = RsSiny. In the 
framework of our approximations (iKopalll 19591) 



b' 



a(l - e ) cos / 
1 + e sin CO 



(16) 



Thus differentiating Eq. ( fT3T l with respect to the time and insert- 
ing equations (flSl l and (fT6] l. we find 



H = 



b^ tani<— >+ 
at 



(2b'^ - R^,) 



ecos w 



<-) 
I + esin CO dt 



(17) 



This equation is valid for all the values of b' except \b'\ - Rs, 
i.e. when the projected trajectory of the planet is tangential to the 
disc of the star and there is no transit. 



2.3. Metiiod to observe tine TDV and its accuracy 

Equation ( fTTb gives the TDV as a function of the variation in the 
orbital elements and parameters derived from the modelling of 
the transit light curve. On the other hand, the variations in the 
orbital elements depend on the orbital parameters of the system, 
the stellar rotation rate, and the inclination of the stellar spin 
to the hne of sight (cf. Eqs. |7] [8] and |9l). The stellar rotation 
period can be derived from the rotational line broadening, the 
radius, and the inclination of the stellar spin i^, or from the ro- 
tational modulation of its optical flux. The inclination can be 
measured in the case of a la te-type star by asteroseismic meth- 
ods (e.g., iBallot et"al]|2006h or, for a moderately active star, by 
a suitable spot modelling applied to high-precision photomet- 
ric time series such as those obtained with CoRoT or Kepler 
(Mosser et al. 2009). 

In the case of a circular orbit, the precession of the line of 
the nodes induces an apparent inclination change that modifies 
the transit duration. The first term on the right hand side of Eq. 
(fTTl i is indeed found to be caused by nodal precession, which 
we denote as [the subscript w indicates that the argument 

of the periastron is fixed during th e variation], is identic al to the 
expressions of dtd/dt found in Miralda-Escude (2002) for cir- 
cular orbits with the advantage that we do not introduce angles 
that cannot be measured because our (i, j) plane is the plane of 
the sky and not the invariable plane of the system. The second 
term in Eq. (fTTl i. which we denote as (^)., is non-vanishing 
only for an eccentric orb it and is related to th e precession of the 
line of the apsides (cf. iPal & KocsisI 12008 ^. As already noted 
by several authors, (^). vanishes when the impact parameter 
b' = y/2RJ2. 

The duration of the transit can be determined with an ac- 
curacy that depends on the depth of the transit, the accuracy of 
the photometry, and the presence of noise sources, e.g., related 
to stellar activity. In the case of the spatial mission CoRoT, the 
photon-limited photometric accuracy ranges typically from ~ 75 
ppm in a hour integration time to ~ 1 130 ppm/hour for a star in 



^ The impact parameter is often given in units of the stellar radius, 
although here we use the non-normalized value. 
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Name 


1 cIT 




e 




CO ( ) 


^ ( ) 


V mag 


Keis. 


/IH (s/yr) 


Art|,^ (s/yr) 


A//|,- (s/yr) 


CoRoT- 1 


5950 


± 150 









77 ± 11 


13.6 


3, 12 


-1.12 


-1.12 




*CoRoT-2 


5625 


± 120 









-7.2 ± 4.5 


12.6 


1, 5 


0.92 


0.92 




T T A '~V Tl '^ 

HAT-P-2 


6290 


± 110 


0.5163 ± 


0.0025 


185.22 ± 0.95 


0.2 ± 12.5 


8.71 


17, 11 


3.04 




3.04 


WASP-33 


7430 


± 100 









-107.7 ±1.6 


8.3 


6 


-34.80 


-34.80 




Kepler- 8 


6213 


± 150 









-26.9 ± 4.6 


13.89 


8, 10 


2.06 


2.06 




TrFS A 




+ 7S 

m I J 


n 

u 






_f. T + zt 7 

— \J.J It / 


1 i .u 


It 


U. / J 


n 7'? 




*WASP-3 


6400 


± 100 









3.3 ± 3.45 


10.5 


16 


-1.03 


-1.03 




WASP- 17 


6550 


± 100 


0.24 ± 


0.07 


278 ± 6.5 


-148.5 ±5.1 


11.6 


2, 15 


-5.82 


-1.16 


-4.66 


WASP- 18 


6400 


± 100 


0.0085 ± 


0.001 


-92.1 ±4.9 


-5.0 ±3.1 


9.3 


7, 15 


1.79 


1.75 


0.04 


*WASP-38 


6150 


± 80 


0.0314 ± 


0.0046 


344.0 ± 17.5 


15 ±43 


9.4 


4, 13 


-0.71 


-0.69 


-0.005 


XO-3 


6429 


± 100 


0.2884 ± 


0.0035 


346.3 ± 1.3 


-37.3 ± 3.7 


9.8 


18,9 


12.70 


12.73 


-0.03 



The asterisk de notes the systems tha t 



Alonso et al. 



Table 1. Transiting planetary systems with measured misalignment for which // > 2 x 10 ^ 

have a significant TDV when taking the extreme possible values for /I, see text for details. References;_J 

2:Rndersonet al. (2010); 3: Barae et al. (2008); 4: Barros et al. (2011); 5: Bouchy et al. (2008); 6: Col Uer Cameron et al. 
7: iHenier et al. (2009); 8: Jenkins et al. (2010b); 9: Johns-Ki-ull et al. (,200 
12 : iPont et al (201(3) ; 13 : Simp son et al. (2010.) ; 14: .Sozzetti et al., (1200' 
n iWinn et al.((l2007h : 18: IWinn e t al. (2009^ 



( 2008 ): 

. (201C): 

10: Kipping & Bakos (20111) ; 1 1: |P al et al.' (20 1(] ); 
15: .Triaud et al., (,201Q) ; 16: iTripathi et al., (,201(31) ; 



the magnitude range 11 < R < 16 dAi grain et al.ll2009h . The 
Kepler telescope can reach ~ 30 ppm/hour to ~ 700 ppm/hour 
over the same range of magnitudes (Jenkins et al.. .2010a ). One 
of the most precise determination of transit duration was that 
obtained for CoRoT-11. In this case, stellar magnetic activity 
has a negligible impact because the stellar optical flux does not 
show a detectable modulation induced by starspots. Therefore, 
the accuracy is dominated by the transit depth and the photo- 
metric accuracy. The depth of the transit is x 0.011 mag and 
the star has V - 12.94 giving a photometric accuracy of ~ 250 
parts per million for measurement s collected b^ the CoRoT tele- 
scope with a cadence of ~ 130 s (iGandolfi et aLlffOK^ . The ac- 
curacy of the transit duration is ~ 50 s considering a sequence 
of ~ 50 consecutive transits for the modelling. Since the aper- 
ture of the CoRoT telescope is 27 cm, a dedicated space-borne 
telescope with an aperture of ~ 3 m should easily reach an ac- 
curacy of ~ 5 s in the measure of H by acquiring data with a 
shorter cadence. Considering a time interval of 10 years, i.e., 
~ 3.2 X 10^ s, we estimate that H ^ 2 x 10"^ can be as- 
sumed as our detection limit, corresponding to a variation of the 
transit duration AH of ~ 0.6 s/yr. Miralda-Escude (2002) ob- 
tained the same limit for HST observations of the brighter target 
HD 209458 (V - 7.64) over a time interval of three years. Hot 
stars (Teff > 6000 - 6200 K) are the most suitable candidates 
for this kind of measurement because of their more likely mis- 
alignment and their low level of magnetic activity, making the 
distortion of the transit profile and its depth variations induced 
by starspots negligible. 

In addition to the light dip produced by the transit of the 
planet across the disc of its host star, it is possible to ob- 
serve the occultation of the planet by the star, i.e., the sec- 
ondary eclipse. This observation is easier in the infrared be- 
cause of the most favourable flux ratio and has indeed been per- 
fo rmed for some of the systems cit ed in this study, e.g ., CoRoT- 
I (lAIonso et al I l2009h . CoRoT-2 (iGillon et al. II20IOI). WASP- 
12 (ICroUetaTTl20Ilh . and WASP- 18 (Nvmever et al.' '20l(k 
The s econdary eclipse of OGLE-TR-56 ( Sing & L6pez-Morales| 
12009 ) has been observed in the visible from the ground. The tim- 
ing of the mid-occultation with an uncertainty of a few seconds 
has also allowed us to estimate e cos co for those systems with a 
precision of the order of 10"^. Moreover, the different durations 
of the transit and the occultation c an be used to measure e si n w, 
although with a lower precision (iRagozzine & Wolfll2009l) . In 
Eq. ( fTTT i. we note that, whereas an increasing inclination of the 



orbit (i.e., di/dt > 0) causes an increase in the duration of both 
the transit and the occultation, the sign of (^). changes in the 
case of the occultation because in this case n must be added to 
the argument of the periastron in Eq. ( fTSl ). 

2.4. Variations in fhe time of the mid-transits for eccentric 
systems 

In the case of an eccentric orbit, the precession of the periastron 
will not only cause TDVs but also a variation in the time of the 
mid-transit (hereafter TTV) leading to a slow change in the ob- 
served orbital period fobs- Since we have neglected the effect of 
the changing inclination of the system / on the true anomaly of 
the planet at the time of mid-transit (cf. Eq. fT4land[T5]). we as- 
sume that the nodal precession will not affect the rate of change 
of the observed period Pobs- Therefore, we use the derivations of 
Pal & Ko csis (2008), valid for / = 90° and arbitrary eccentricity 



(1 



2)3/2 



p2 

e cos CO r. 



ob.s 



1^^) 

(1 -H esinw)^ n \ dt 



(18) 



where Pq is the true orbital period, i.e., Pq = 2nln. The diff'er- 
ence O - C between the observed epoch of mid transit and that 
computed with a constant-period ephemerides is 



0-C 



1 



(19) 



where is the number of transits elapsed from the initial 
epoch and we have considered a time interval much shorter than 
2nl{^) so that Pobs, given by Eq. (fTSl ), is virtually constant. 
We note that the sign of Pobs is the opposite for the occultation 
because in this case n must be added to the argument of the peri- 
astron u) in Eq. ( fTsT l. Therefore, the (9 - C of the occultation will 
have the opposite sign of that of the transit, by analogy with the 
case of eclipsing clo se binaries with eccentric orbits (cf., e.g., 
iGimenez et 571119871) . 

Finally, we note that the presence of a third body in the sys- 
tem can induce a further precessio n of the orbital plane leading 
to an additional term in (cf. iMiralda-Escude 2002) . that 
will change the value of H and Fobs- If the semi-major axis of the 
third body orbit is much larger than that of the planet, the light- 
time effect due to the motion of the star-planet system around the 
barycentre of the triple system is the dominant eff'ect and we ex- 
pect to observe the same variation in the times of both the transit 
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and the occultation. If the third body orbit is close, it is Hkely to 
be in a mean motion resonance with the transiting planet orbit, 
a configuration which ensures the stability of the system. In this 
case, TTVs of several tens of seconds or minutes are expected, 
even for a th ird body with a mass comparable to that of the Earth 
dAgol et al. I I2OO5'). Moreover, radial velocity measurements can 
be used to look for a close third body, allowing us to test this 
interpretation. 

In conclusion, the observation of a TDV without an accom- 
panying TTV can be considered as an indication that nodal pre- 
cession is the dominant effect, at least for a circular orbit. In the 
case of an eccentric orbit, an additional contribution to the TTV 
arises from the precession of the periastron. In principle, it can 
be distinguished because it produces opposite TTVs for the tran- 
sit and the occultation. 



Name 



r,ff (K) V mag Refs. 



CoRoT-11 


6440 ± 


120 


12.8 


GalO 


-48.05 


HAT-P-6 


6570 ± 


80 


10.4 


N0O8 


-1.04 


HAT-P-9 


6350 ± 


150 


12.3 


Sh09 


-1.17 


OGLE-TR-56 


6119± 


62 


16.6 


T0O8 


-0.67 


0GLE-TR-L9 


6933 ± 


60 


13.97/ 


Sn09 


-38.59 


WASP- 19 


5500 ± 


100 


12.3 


Hell 


-2.59 



Table 2. Transiting planetary systems with undetected misalign- 
ment and eccentricity with H > 2 x 10"^ for e < 25°. The 
value of AH is given in seconds per year an d is computed 
for U = 90° and e = 25° . Refe rences: GalO: iGandolfi et al.1 
('2010'); Hell:'Hellieretal.' ('2011'); La09 'Latham et al.' ("20091); 
N0O8: Novesetal. (2008); Sh09: Shporer et aL (2009); Sn09: 
ISnellen et al.l (12009h : T0O8: ITorres et al.l (120081) . 



3. Applications 

3.1. Systems with measured misalignment 

At the time of writing, the projections of the obUquity of 3 1 tran- 
siting systems with hot Jupiters have been measured. However, 
this measure is usually subject to an uncertainty of 50-100 per- 
cent when \A\ < 10°. Assuming H - 2x 10~^ as the detection 
limit, there are eight known misaligned systems for which the 
TDV can be measurable. In Table[T] we list their relevant param- 
eters, from the left to the right: the name of the star, its effective 
temperature, the eccentricity of the orbit, the argument of peri- 
astron, the V magnitude of the star, the references, the change 
of the duration of the transit in seconds per year A//, the varia- 
tion in the duration caused by nodal precession A//|„ in seconds 
per year, and, for eccentric systems, the variation in the dura- 
tion caused by apsidal precession A//|,, assuming /j = 90°, also 
in seconds per year. We note that there are three more systems, 
i.e., CoRoT-2, WASP-3, and WASP-38, denoted by an asterisk 
in Table [T] for which the TDVs can be measured if we assume 
the extreme possible values for their measured A. We note that 
for HAT-P-2, / - 90° so that b' - Q within the uncertainties and 
the nodal precession term vanishes at the present time. 

A significant departure from the predicted TDV value can be 
inte rpreted as an indication of either a third body in the system 
(cf. lMiralda-Escu"dgl2002h or an oblique planetary spin axis, if 
the system is very young. Specifically, if the tidal quality factor 
of the planet Q'^ ~ 10^ - 10^, as predicted by some models for 
coreless massive planets (cf.. e.g.. iGoodman & Lacknerll2009t 
IOgilvidl2009l ; |Papai oizou & Ivanov 1120 1 Oh . the tidal damping of 
the planetary obliquity takes place over timescales of the order of 
10 - 100 Myr that would lead to an observable additional nodal 
precession. This happens because an additional term is present 
in Eq. (|7]l when the planetary spin is not aligned with the or- 
bital angular momentum, with a coefficient Ap ~ 100 A, that 
gives rise to a significant contribution to the variation in the in- 
clination of the orbital plane. The detection of this effect would 
provide a lower limit to the largely unkno wn planetary tidal fac- 
tors in the case of young sys tems (cf., e.g., lOgilvie & Linll20()4t 
IGoodman & Lacknerll2009l) . 

3.2. Systems witii a circular orbit 

To date, 32 transiting exoplanetary systems with undetected ec- 
centricity and misalignment are known. If the eccentricity of the 
orbit is close to zero, then the only significant source of variation 
of the duration of the transit is the spin-orbit misalignment. We 
give in Fig.[T]the number of systems for with \H\ > 2 x 10"^, 



i.e., expected to be detectable within three to ten years of obser- 
vations, computed on a grid of 5° resolution for 0° < e < 85° 
and 0° < is < 90°. The most favourable case occurs for the 
greatest value of the obliquity and small values of i^, because the 
angular velocity of the star is higher than the projected measured 
value, at least for inclination values that do not lead to a rotation 
rate exceeding the break-up velocity. The other extreme case is 
more interesting, as it shows that small obliquities can indeed 
be constrained by measuring TDVs. We list in Table|2] from left 
to right, the name, the eff'ective temperature, the V magnitude, 
the reference for the orbit, and AH in seconds per year of the 
six stars for which a value e < 25° can be detected assuming a 
2 X 10"^ level of accuracy for H. The three systems that have the 
fastest nodal precession, even for a quasi-aligned configuration, 
are CoRoT-11, OGLE-TR-9, and WASP-19 for which e < 10° 
can be measured provided that 70° < < 90°. 

In the case of stars hotter than T^ff ~ 6500 K, an accu- 
rate measurement of A through the Rossiter-McLaughlin ef- 
fect can be difficult owing to their rapid rotation and relative 
pau city of spectral lines (cf., e.g., the case of CoRoT-11 in 
Gan dolfi et a In those cases, an upper limit to e can be 

derived solely from the measurement of the nodal precession be- 
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Fig. 1. The current number of systems with a circular orbit ex- 
pected to have H >2x 10"^ as a function of the obliquity e and 
the inclination of the stellar spin to the line of sight i^, in bins of 
5° x5°. 
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cause sin 4 sin /I < 1, so that Eq. (|7]i gives 



cos e > 



AM} 



<-) 



(20) 



where the variation in the incHnation is derived from the TDV 
and the parameters of the transit (cf. Eq.fTTb. 

3.3. Systems with eccentric orbit 

In the case of eccentric systems, the duration of the transit will 
vary, even if the system is aligned, owing to the apsidal pre- 
cession. Out of 28 eccentric systems for which no significant 
Rossiter-McLaughlin effect has been reported, there are three 
that reach \H\ > 2 x 10"^. Table |3] lists their names, the effec- 
tive temperature of the star, the orbit eccentricity, the argument 
of the periastron, the V magnitude, the reference, the periods 
of apsidal precession caused by the relativistic effect Pgr^ the 
tidal bulges raised on the star and the planet ftid, the oblate- 
ness of the star and the planet induced by their rotation Pj^, and 
AH for is = 90° in seconds per year. For all of them, the pre- 
dominant part of the apsidal precession is caused by the tidal 
bulges, and for HAT-P-23 and WASP- 12, the oblateness effect is 
st ronger than the re l ativist ic effect, in contrast to the assumptions 
of lJordan & BakosI (l2008l) . However, this does not hold true for 
aU the eccentric systems and in general we find that Pgr ^ Pj, 
and Pgr ~ Ptid- We note that for small obliquities, the rate of 



Name 



HAT-P-23 



HAT-P-24 



WASP-12 



Teff (K) 


5905 ± 80 


6373 + 80 


6250 ± 150 


e 


0.106 + 0.044 


0.067 ± 0.024 


0.049 + 0.015 


ojC) 


118 ±25 


197 ± 36 


286 + 0.11 


V mag 


11.94 


11.82 


11.69 


Refs. 


BalO 


BalO 


He09, Hull 


^-GR (yr) 


2272 


12047 


1705 


^-iid (yr) 


133 


6407 


18 


P,. (yr) 


1717 


50243 


272 


AH (s/yr) 


15.10 


1.26 


-13.42 



Table 3. Eccentric transiting planetary systems with undetected 
misalignment with \H\ > 2 X 10~^. The value of AH is com- 
puted for = 90°. Refere nces: BalO:lBakos et all (l2010h : He09: 
iHebb et all dlool : Hul l: lHusnoo et all (1201 Ih . 



change in the transit duration does not vary significantly, and 
only two additional systems reach the level of \H\ > 2 x 10"** 
for 10° < e < 25°, namely HAT-P-14/WASP-27 and HAT-P-21. 
We conclude that the measurement of TDVs to detect obliquity 
is the most interesting for systems with a circular orbit. 



4. Systems with observed TTV 

The TTVs of some transiting systems have already been de- 
tected, although in some cases with an accuracy lower than con- 
sidered above. The cause of the effect can be the presence of a 
third body in the system (e.g., Miralda-Escude 2002; Agol et al. 
l2005h . a chan ge in the orbital pe riod owing to the tidal interac- 
tion (cf., e.g.. iLanza et aljr201 ll) or, in the case of an eccentric 
orbit, the apsidal precession. 

Instances of TTVs have been detected for the system WASP- 
3 (^acieiewski et al. 2010). Even considering the highest al- 
lowed eccentricity, apsidal precession alone cannot explain the 
measured value. Furthermore, the period found in the O - C's 



is three orders of magnitude shorter than the period of the pre- 
cession, suggesting that the TTVs are due to a third body in the 
system. 

Finally, we consider the system WASP-12, for which apsi- 
dal precession, in contrast to the other one, could account for 
the TTVs. The eccentricity of this system is quite uncertain, 
and was reported to be e - 0.049 ± 0.015 in the discovery 
paper of lHebbet all (l2009l) . but later constrained to the much 
smaller value e = 0.017+°o!i (iHusnoo etalJl201lh . Two tran- 
sits, which are 661 and 683 epochs away f rom the initial epoch 
of the ephemerides of Hebb et al.l (|20QSj), were observed by 
Macieie wski et a n(l201lV with a formal photometric error in the 
range between 0.6 and 0.7 mmag for the individual measure- 
ments. Assuming the discovery paper value for the eccentricity, 
we estimate an O - C ~ 140 s for their epoch 683, which is 
incomp atible with the observ ed value. Adopting the eccentric- 
ity after iHusnoo et alJ (l201 ll) . we find O - C ~ 42 s, in excel- 
lent agreement with their measured value for their epoch 683. 
However, the predicted O - C disagrees with that observed at 
their epoch 661. Therefore, additional high-precision photome- 
try of this system would be required to reach a definite conclu- 
sion. 



5. Conclusions 

We have shown that the variation in the transit duration should 
be readily measurable for several known planetary systems, and 
can be used to constrain the value of the angle between the spin 
axis of the star and the orbital angular momentum. Assuming 
that the planet is in a quasi-synchronous and aligned state with 
respect to its orbit, and neglecting tidal interaction and a possi- 
ble third body, we have derived the rate of change in the orbital 
elements produced by the deformation of the bodies resulting 
from the centrifugal flattening and the tidal bulges as well as 
the general relat i vity pr eces sion. After combining the results of 
lMira.lda-Escudi (l2002h and iPal & KocsisI (l2008l) . we obtained 
a general expression for the rate of change in the transit du- 
ration, which is valid for eccentric and/or misaligned systems. 
Assuming an inclination of the stellar spin axis = 90°, a rela- 
tive photometric accuracy of about 100 ppm in a hour integration 
time, a typical transit depth of 0.01 mag, and an observation in- 
terval of 3 - 10 years, we found that: 

- eleven systems should display observable TDVs if their obli- 
quity is consistent with the values of A measured through the 
Rossiter-McLaughlin effect; they are: CoRoT-1, CoRoT-2, 
HD147506/HAT-P-2, HD15082/WASP-33, Kepler-8, TrES- 
4, WASP-3, WASP-17, WASP-18, WASP-38, and XO-3. 

- six systems should display observable TDVs if their ec- 
centricity is close to zero and their obliquity smaller than 
25°; they are: CoRoT-1 1, HAT-P-6, HAT-P-9, OGLE-TR-56, 
OGLE-TR-L9, and WASP-19. 

- three systems may display observable TDVs, given the un- 
certainty in their measured eccentricity; they are: HAT-P-23, 
HAT-P-24, and WASP-12; moreover, HAT-P-14/WASP-27 
and HAT-P-21 could have observable TDVs, if their obli- 
quity is 10° < e < 25°. 

For circular systems, an upper limit to the obliquity can be 
obtained using TDVs, which can be particularly interesting in 
the case of hot host stars whose rapid rotation gives rise to a size- 
able quadrupole potential and rapid nodal precession. In princi- 
ple, tidal or gravitational interaction with another body and nodal 
precession effects can be separated because the former two pro- 
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duce TTVs and TDVs, whereas the latter gives rise to TDVs only 
(cf. Sect. |23]for details). 

The TDVs can be used to constrain the obliquity of the 
star (and the planet in systems younger than 10 - 100 Myr), 
and the J2 of the star and the planet, so that they can give in- 
formation about the formation scenarios, or constrain the star 
and planet's internal structure through the apsidal motion con- 
stant or Love number Moreover, they can provide an indica- 
tion of a third body, even of a mass comparable to that of the 
Earth, because this body would make an additional contribu- 
tion to the nod al pre cession or TTV ( see Sects. 12.41 14] and 
iMiralda-Escude 2002': 'Ago l et aL 1 12005). As of now, the spa- 
tial mission Kepler yields the highest cumulative precision and 
TTVs h ave been used to con firm the presence of multiple plan- 
ets (e.g. lHolman et alj|2010h . However, the determination of the 
transit durations of the non-multiple systems does not reach the 
level of accuracy discusse d here. In their independent analysis 
of Kepler-4b to Kepler-Sb. lKipping & Bakos (.201 1.) obtained an 
accuracy of several tens of seconds on the possible TDVs and 
TTVs for those planets using only the discovery photometry, 
which spans about 40 days. However, Kepler has been in op- 
eration since May 2009, meaning that the number of observed 
transits should have increased by a factor of about 40 at the end 
of the mission, leading to a significant improvement in the accu- 
racy of the measured transit durations. In the case of Kepler-4b, a 
2 cr marginal detection of eccentricity has been claimed, but only 
an upper limit at e < 0.43 can be assumed to 95% confidence. 
U sing the parameters inferred from the eccentric fit to Kepler-4b 
of lKipping & BakosI ([201 1), the apsidal precession would cause 
a TDV of only about 0. 17;^^^^ s over the two years of the mission 
for e - 0.25, owing to the unfavourable orientation of the orbit 
with respect to the line of sight (co - 84.5;^[g j). Another inter- 
esting system in the Kepler field of view is HAT-P-7. However, 
even assuming the most favoura ble value of the obl iquity com- 
patible with the observations by I Winn et all (l2009bl) . we do not 
expect a TDV variatio n greater than one se cond over a time span 
of ~ 15 yr Moreover. IWinn et al.l (l2009bl) present evidence of a 
third body in the system, thus a detection of a TDV by Kepler 
might confirm that a companion exists. 

All other single planets discovered by Kepler so far have pos- 
sible TDVs ranging from a few tenths of to one second over two 
years of observation, according to their different parameters. By 
the end of the mission, the accumulated photometry will surely 
be sufficient to place stringent constraints on the models param- 
eters, and if the systems are misaligned, allow us to probe their 
stellar interiors through their TDVs. 
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